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We propose a measurement protocol to generate path-entangled NOON states conditionally from 
two pulsed type II optical parametric oscillators. We calculate the fidelity of the produced states 
and the success probability of the protocol. The trigger detectors are assumed to have finite dead 
time, and for short pulse trigger fields they are modeled as on/off detectors with finite efficiency. 
Continuous-wave operation of the parametric oscillators is also considered. 
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I. INTRODUCTION 

Nonclassical states of light have many applications, 
and a number of different protocols exist for the genera- 
tion of various classes of states. The two-mode maximally 
entangled A^-photon states 



|NOON) = ^ (|7V, 0) + e*'^|0, N)) 
v2 



(1) 



the so-called NOON states, are particularly interesting 
because a single-photon phase shift of x induced in one of 
the two components changes the relative phase of the two 
terms by Nx- This special property of NOON states may 
be utilized to enhance spatial resolution in (quantum) mi- 
croscopy and lithography [ij , and in interferometry it has 
been shown that a certain measurement strategy, using 
NOON states, leads to a phase estimation error scaling 
as L~^^'^N^^^^ if the phase to be estimated is known to 
lie within an interval from —it/ L to tt/L, where Nt is the 
total number of photons used in the measurements 0. 

— 1/2 

This is better than the classical shot noise limit of Nj, ' , 
and NOON states are thus useful to perform accurate 
measurements and may be a valuable field resource in 
sensors. NOON states are also a source of entanglement 
with applications in quantum information protocols and 
in fundamental studies such as tests of Bells inequality 
i- 

It is thus of great interest to be able to produce NOON 
states, and various NOON state generation schemes have 
been suggested theoretically 0,^5, 6, 7, 8, 9] and studied 
in experiments 0, i, [M [Hi- The TV = 1 and = 2 
NOON states may be generated by combining either a 
single photon and a vacuum state or two single-photon 
states on a 50 : 50 beam splitter, but this simple approach 
is not directly extendable to > 2, and we shall thus 
mainly be concerned with generation of = 3 NOON 
states in the present paper, even though the suggested 
protocol is, in principle, applicable for all A^. Mitchell, 
Lundeen, and Steinberg have generated NOON states 
with A^ = 3 from a pair of down converted photons 
and a local oscillator photon using certain pola rization 
transforming components and post-selection • In this 
experiment, however, the successful generation of the 



NOON state is witnessed by a destructive detection of 
the state. In the present paper we propose and analyze 
in detail a nondestructive generation protocol, which con- 
ditions the successful generation of the A^-photon NOON 
state on the registration of A' photo detection events in 
other field modes, and which uses as resource only linear 
optics and the output from two optical parametric oscil- 
lators (OPOs). The protocol does not rely on efficient 
photo detection. The analysis is carried out in terms of 
Wigner function formalism, and effects of finite detector 
efficiency and finite detector dead time are considered. 

Conditional generation of nonclassical states occupying 
a single mode has been investigated both experimentally 
and theoretically [12, IS, ii, 15, 16, 17, 18, 19, 20]. With 
the correlated output from a single nondegenerate OPO 
it is, for instance, possible to generate n-photon Pock 
states of light in the signal beam conditioned on n photo 
detections in the idler (trigger) beam [H, [H, and 
in principle the entanglement of a highly squeezed two- 
mode field from an OPO makes it possible to prepare any 
state in the signal beam that can be either measured as 
an eigenstate of a suitable observable of the idler beam or 
produced as the final state of a generalized measurement. 
The basic idea of the protocol proposed in the present 
paper is to mix the output from two OPOs and employ 
the entanglement to prepare a two-mode state in two of 
the output beams by detection of the desired output state 
in the remaining beams. 

In Sec. II we explain the NOON state generation pro- 
tocol in detail. In Sec. Ill we analyze the performance 
of the protocol quantitatively for pulsed OPO sources. 
We provide the fidelity of the generated states and the 
success probability. In Sec. IV we consider production of 
NOON states from continuous-wave OPO sources, and 
Sec. V concludes the paper. 



II. EXPERIMENTAL SETUP FOR NOON 
STATE GENERATION 

The experimental setup is illustrated in Fig. [TJ Two 
pulses of two-mode squeezed states are generated by two 
identical OPOs via type II parametric down conversion. 
The field mode operators of the modes generated by the 
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first OPO are denoted a+ and a_, respectively, while the 
field mode operators of the modes generated by the sec- 
ond OPO are denoted 6+ and 5_, respectively. For def- 
initeness, we assume that the plus modes are vertically 
polarized and that the minus modes are horizontally po- 
larized. The modes are separated spatially by the first 
two polarizing beam splitters, and the third polarizing 
beam splitter combines the a_ and modes, which are 
subsequently subjected to the NOON state measurement 
proposed in and illustrated for iV = 3 in Fig. [TJ The 
idea behind this measurement is to apply the highly non- 



linear operator An = a_ — ) to the state. The 

result is only nonzero if either the a_ mode or the 6+ 
mode contains at least N photons. On the other hand, 
if the squeezing is sufficiently small, it is unlikely to have 
more than a total of N photons in the two trigger modes, 
and by conditioning on the successful application of An, 
we select the pulses of the system where TV photon pairs 
are generated in one OPO and zero photon pairs in the 
other. It is equally probable that the photons originate 
from the first OPO or from the second OPO, and, as we 
shall see in detail below, the result is that a NOON state 
is generated conditionally in the output modes a_|_ and 

As stated in [2l[, An can be rewritten as a simple 
product of single photon annihilation operators 



N ^ 



i0\ i2Tm/N 



(2) 



and it is thus possible to implement An hy means of 
beam splitters and photo detectors. We first consider odd 
values of N. Beam splitters are used to divide the input 
into N distinct spatial modes labeled by n = 1, . . . , TV. 
The beam splitter reflectivities are chosen to obtain the 
same expectation value of the intensity in each of the 
modes. The vertically polarized modes are then phase 
shifted by the factor g^^nn/N+nr relative to the horizon- 
tally polarized modes, i.e., 6+ — fe+e'^'^"/^, and finally 
polarizing beam splitters with principal planes oriented 
at 45° relative to the horizontal polarization transform 
a_ and -fo+e'^^/jv j^^^ (a_-S+e*2™/^)/V2 (the trans- 
mitted mode) and (a_ + 6+e'2^"/^)/V2 (the reflected 
mode) [l^l. The annihilation of a photon in each of 
the modes transmitted by the beam splitters witnesses 
the overall application of the operator An- If one ob- 
serves both reflected and transmitted modes simultane- 
ously, one conditions on detection events in all the trans- 
mitted modes and no detection events in all the reflected 
modes. If detection events are instead observed in all the 
reflected modes and in none of the transmitted modes, 
an operator of the form ^ is also obtained, but 9 is ef- 
fectively transformed into 9 + tt due to the phase shift 
at the polarizing beam splitter, and the value of (f> of the 
generated NOON states is changed by Ntt (see below). 
The success probability is thus increased by a factor of 
two if both outcomes are accepted. 




FIG. 1: Experimental setup for NOON state generation. 
OPO, optical parametric oscillator; PBS, polarizing beam 
splitter; PS, phase shifter; and APD, avalanche photo diode. 
The part of the setup enclosed in the dashed box performs the 
NOON state measurement, and here it is shown for A'^ = 3. 
Note that the polarizing beam splitters inside the box are ori- 
ented at 45°. The numbers denote beam splitter reflectivities 
of 1/3 and 1/2, and the three phase shifters transform b+ 
into —b+e^'^^"^^, where n = 1,2,3, respectively. See text for 
details. 



For even values of iV a similar measurement scheme 
is applicable, but it is sufficient to divide the field into 
N/2 spatial modes initially, and in this case the NOON 
state generation is conditioned on detection events in 
both transmitted and reflected modes (see [2l|). 



III. PERFORMANCE OF THE PROTOCOL 

After this presentation of the basic idea and the phys- 
ical setup we now consider the actual outcome of the 
detection process. For short pulse OPO output the dead 
time of the photo detectors may typically be longer than 
the pulse duration, and we shall thus assume that it is 
impossible to obtain more than a single detection event 
per detector per pulse, i.e., if the detector efficiency is 
unity, the detectors are only able to distinguish between 
vacuum and states different from the vacuum state. Such 
detectors are denoted ori/off detectors, and they are dis- 
cussed in detail in Ref. [231 . The finite dead time of the 
detectors is not severe to the measurement procedure de- 
scribed in [m because the on/off detector model and the 
conventional photo detector model, represented by the 
annihilation operator, lead to identical signal states if the 
total number of photons in the idler modes is guaranteed 
to be less than or equal to the number of conditioning 
detection events, i.e., N. 

We analyze the performance of the setup using Gaus- 
sian Wigner function formalism [l^, [l^, |20| | , which is 
applicable because the squeezed states generated by the 
OPOs and the vacuum states coupled into the system 
via the beam splitters are all Gaussian. In general, the 
Wigner function of an n-mode Gaussian state with zero 
mean field amplitude takes the form 



Wvixi,pi, . . .,Xn,Pn) 



1 



(3) 
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where y = {xi,pi, . . . ,Xn,Pn)^ and V is the 2n x 2n 
covariance matrix. If Ci denotes the field mode annihi- 
lation operator of mode i, the elements of V are given 
in terms of the real and imaginary parts of the expecta- 
tion values (cjcj) and (ciCj). Note that for a multi-mode 
Gaussian state we are free to include only the modes of 
interest in ^ because the partial trace operation over un- 
observed modes is equivalent to integration over the cor- 
responding quadrature variables. A unit efficiency 'on' 
detection in mode i projects mode i on the subspace of 
Hilbert space that is orthogonal to the vacuum state, i.e., 
the Wigner function is multiplied by (1 — 2TTWo{xi,Pi)), 
where Wo{x,p) = exp(— a;^— p^)/7r is the Wigner function 
of the vacuum state, the variables Xi and pi are integrated 
out, and the state is renormalized. Since the Gaussian 
nature of a state is preserved under linear transforma- 
tions, and since a detector with single-photon efficiency 
ry is equivalent to a beam splitter with transmission rj fol- 
lowed by a unit efficiency detector [2^ , effects of non-unit 
detector efficiency are easily included in the covariance 
matrix. 



To calculate 



and 



explicitly we note that 



the state generated by the OPOs is [2j 



IV'i> = (l-r^)EE 



r"+™|n,n, m, m) 



(4) 



n— m—0 



where r is the squeezing parameter and the modes are 
listed in the order: a+, a_, &+, We assume that N is 
odd and consider the transmitted trigger modes (which 
we number from 1 to N), the a+ mode (mode N + 1), 
and the 6_ mode (mode iV + 2). By expressing the field 
operators of the trigger modes (those observed by the 
unit efficiency detectors) in terms of a_, 6+, and field 
operators representing vacuum states we find 



^27ri(k-j)/N 



27V' 



(5) 



where j e {1, 2, . . . , N}, G {1, 2, . . . , iV}, A = wVi^ " 
r^), and we allow of a constant phase shift 9 of 6+ relative 



to a_. Furthermore 

(cjv+iCAT+i) = {ijM\a'\_a+\ip^) = rV(l - r^), (6) 
(4+25^+2> ^{HblU^^)=ry{l-r^), (7) 



7] r 



(8) 



2iV 1 - H 



^2Trik/N+ie 



(9) 



and 



[CjCk) = [CN+lCN+l) = {CN+2CN+2) = {CN+1CN+2) = 

{cN+iCN+2) = (4cw+i> = (4cw+2> = 0. (10) 

For even values of N the factors yJrj/{2N) are replaced by 
y^rj/N. Note that loss in the signal beam may be taken 
into account by performing the transformations a+ —> 
yjrjlaj^ and 6_ — > ,Jrjlb^ in the above expressions, where 
1 — is the loss. 

The NOON state fidelity F/v of the signal state con- 
ditioned on N photo detection events in the transmitted 
trigger modes is 

Fn = —— I W-t^oOT>i{xN+l,PN+l,XN+2,PN+2) 
J 

(j[il-^^Woix,,p,))^ 

/N+2 \ 
Wv{xi,Pl, . . . ,Xn+2,PN+2) Y\. d^i^Pi ' (11) 

where IVnoon is the Wigner function of the NOON state 
dH), and 



Pn = J (jl{l-27rWQ{x,,p,)) 



Wv{xi,Pl, . . . ,Xn+2,PN+2) 



/N+2 \ 

Y[ dx.dp, , (12) 



is the success probability, i.e., the probability to obtain 
the conditioning detection events and produce the NOON 
state in a given pulse of the OPO system. We expand 
the product 



N 



N 



n (1 - 2^M^o(x„pO) = E X{{-'^^Wa{x^.P^)f\ 
i— 1 d i—1 

(13) 

where the sum is over all d = (di, ^2, ■ • ■ , c^tv) with 
di £ {0, 1}, and define the diagonal matrix Jd = 



4 



diag((ii, di, ^2, ^2, ■ • ■ , ^at, dA?) and the nx n identity ma- 
trix /„. Furthermore, we divide the covariance matrix 
into four parts 



V = 



(14) 



where Vu is the 2N x 27V covariance matrix of the trigger 
modes, Yss is the 4x4 covariance matrix of the signal 
modes, while Vts contains the correlations between the 
trigger and the signal modes, and we define the vector 
ys = {xn+i,Pn+i,xn+2,Pn+2V and the matrix 



(15) 



This allows us to write Eqs. p4|) and in the following 
compact forms [1^ 



Fn — 



1^ V 



\/det{l2N + JdVtt) 



and 



Pn = 



E 



(-2)5:f^i'^' 



v/det(/2Ar + JdVtt) 



(16) 



(17) 



Since Wnoon is a product of a polynomial and a Gaus- 
sian the integral in Eq. (jl6p may be evaluated analytically 
and for N = 3 and ?/ = 1 we find 



^3 



(1 - r^y{2 - r^y{3 ~ 2r^)(6 - 5r^) 
18(4- 3r2) ^ 



(18) 



where the optimal value cf) ~ N9 + tt + 2iin, n ^ Z , 
is assumed. Expressions for P/v are given in table [J for 
TV = 1, 2, 3, and 4, and F/v and Pa? are plotted for = 3 
in Figs. [Hand [31 respectively. We observe that high prob- 
abilities are only found in the parameter regime, where 
the fidelity is low. If, for instance, we want a NOON 
state fidelity of at least 0.9, we choose r — 0.14, and if 
77 — 0.25, P3 is of order 10~*. With a repetition rate 
of order 10® s~^ (see [l^) one state is produced every 
second minute on average. The production rate is very 
dependent on detector efficiency, and if rj is increased to 
unity, the rate is increased by approximately a factor of 
60. 

For odd values of N we may observe both reflected 
and transmitted trigger modes and condition on detec- 
tion events in all the transmitted trigger modes and no 
detection events in all the reflected trigger modes, or, vice 
versa. In this case we also include the reflected trigger 
modes in the covariance matrix, which we now denote by 
V'^ . By a similar analysis as above we obtain the success 
probability 



2E 



2^'(-2)2:f=irf. 



d \lAet{hN + J+V+) 



(19) 



iV 


Pn 


i^rx 


1 


\ 

A+l 


■2\ 
(A+1)^ 


2 


(A+1)^ 




3 


A-'(A+4) 


2A-'(3A+4) 


(A+2)^(A+3)(A+6) 


(A+l)^(A+2)^(2A+3)(5A+6) 


4 


A*(A^+6A+6) 




(A+1)2(A+2)^(A^+8A+S) 



TABLE L Success probabilities calculated from Eqs. ((17} and 
dUl). \ = rir'^/{l~r'^). 




FIG. 2: NOON state fidelity F3 (solid lines) and F3+ (dashed 
lines) as a function of squeezing parameter r for r] = 1 (upper 
lines), rj = 0.25 (middle lines), and 77 — > (lower lines). Note 
that in the latter case F3 — F.t . 



where 



diag(di, di, ... ,dN, (In, 1,1, ... , 1, 1), while 
the NOON state fidelity is given by Eq. (HH) with V 
replaced by the matrix 



V^iV+fiV+j, + l2N)-% 



R ' 



(20) 



where is the covariance matrix of the reflected trig- 
ger modes and consists of the correlations between 
the reflected trigger modes and the signal and transmit- 
ted trigger modes. Explicit results for P^ are given in 
table [J for iV = 1 and 3. F^ and P-^ are compared to 
P3 and P3 in Figs. [2] and[3l and it is observed that F^ 
and P^ are both larger than P3 and P3 if r is not large 
(and 77 > 0). For r ^ 1, because in this limit it 

is very unlikely to obtain no detection events in all the 
reflected or in all the transmitted trigger modes. 

In the limit of very small detector efflciency a simple 
expression for the NOON state fidelity for the case of N 
trigger detectors is easily derived without using Wigner 
function formalism. In general, if the state of interest 
is expressed in the photon number basis, the mathemat- 
ical operation corresponding to an 'on' detection is to 
multiply each term by -y/l — (1 — 77)", where n is the 
number of photons in the mode observed by the non- 
unit efficiency detector, trace out the detected mode, 
and renormalize. If 77.77 ^ 1 ^ot all contributing terms. 



(1-77)' 



777/ CX 



and the on/off detector 
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FIG. 3: Success probability Ps (solid lines) and Pj^ (dashed 
lines) as a function of squeezing for rj = 1 (upper lines) and 
■q — 0.25 (lower lines). The dotted lines represent the approx- 
imate expression (|23|) . 



state fidelity is Fi^s = (1 — r^)^, and the success probabil- 
ity is Pi^s — A/(A + 1). Choosing squeezing parameters 
such that Fi^s = Fi, we find that = (4/3)Pi^s in the 
high fidelity limit. It is thus possible to achieve a higher 
success probability using the scheme with two OPOs, but 
the price to pay is a more technically involved setup, 
and NOON states with two different values of are pro- 
duced. For N = 2 the present protocol and combination 
of two single-photon states on a 50 : 50 beam splitter, 
each produced conditionally from a single OPO, lead to 
identical fidelities and success probabilities. Finally we 
note that the photo detector model underestimates Fn 
for 7y > because 1 — (1 — 77)" = 'r]J2"~Q{l — rjY < nrj 
for n = 2, 3, . . . while 1 — (1 — rj)" = nq for n = 0, 1, i.e., 
the 'wrong' terms containing more than N photons are 
given a too large weight. This is also what we observe in 

Fig.m 

In the limit of small r and for odd values of N the 
success probability is given approximately by the simple 
expression 



model becomes equivalent to the photo detector model. 
In this case the density operator of the output state is 
obtained as 

00 00 

p=0 ij=0 

_ (l-r2)^+2 

^ oc 00 I 
^n=Nm=0 ^ '' 

(21) 

where Af is a normalization constant and the traces are 
over the a_ and 6+ modes. This leads to the NOON state 
fidelity 

= (NOON|p|NOON) = (1 - r'^)^+'^, (22) 

where again (j) = NO + tt + 27rn, n G Z, is assumed. 
It is interesting to compare this result with the fidelity 
(1 — 7.2)JV-(-i obtained for production of A^-photon states 
from a single two-mode squeezed state by conditioning 
on A'^ detection events in the idler beam and using detec- 
tors with very small efficiency. If a single-photon state is 
produced by this method and transformed into an = 1 
NOON state as explained in the Introduction, the NOON 



{a^_ - ihe'Tm) = (^^"^ odd). (23) 

Again r]/{2N) must be replaced by rj/N to obtain Pn for 
even values of N. The approximation to P3 is shown in 
Fig.IS 



IV. NOON STATES FROM 
CONTINUOUS- WAVE OPO SOURCES 

Our protocol is not limited to pulsed fields, and for 
completeness we now briefly consider NOON state gen- 
eration from continuously driven OPOs. We assume 
N = 3. For continuous-wave fields each of the three 
detected trigger beams and the two signal beams are 
described by time dependent field operators Ci{t). The 
trigger detections take place in particular modes local- 
ized around the three detection times td, tc2, and fc3, 
and we want to determine the NOON state fidelity of 
an output state occupying one temporal mode in each 
signal beam. Following the general multimode formal- 
ism in Refs. [l^, [2OI, the five relevant modes are spec- 
ified by the mode functions fi{t), and the correspond- 
ing five single mode operators are c; = J fi{t)ci{t)dt. 
In general, we are free to choose the two output mode 
functions at will, and in the present case it is natural to 
choose the mode function which gives rise to the largest 
three-photon state fidelity when three-photon states are 
generated from a single type II continuous- wave OPO. 
Since we are mainly interested in the parameter region 
where the squeezing is small and the NOON state fidelity 
is large, we use the optimal three-photon state mode 
fimction derived for very low beam intensity in (20| . i.e., 
hit) = fbit) = ELi Ck\/jj2exp{-j\t - tck\/2), where 
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FIG. 4: NOON state fidelity as a function of e/7 for states 
generated from a pair of continuous-wave OPO sources when 
conditioning on tiiree simultaneous trigger detection events 
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FIG. 5: Fidelity of NOON states from continuous-wave OPO 
sources as a function of separation between trigger detection 
events {tcs — tci)j for A'' = 3, tcs - ic2 = tc2 - id, and e/7 = 
0.01. 

the coefficients Ck are functions of the intervals between 
the detection times and 7 is the leakage rate of the OPO 
output mirror. We furthermore assume that the trigger 
mode functions are nonzero only in an infinitesimal time 
interval centered at the detection time, which is valid if 
the trigger detections take place on a time scale much 
shorter than 7"^. Since we consider a low intensity con- 
tinuous beam, and since we formally assume that the 
trigger detectors only register the light field in infinitesi- 
mal time intervals around the detection times, the anni- 
hilation operator detector model is perfectly valid in this 
case and detector dead time is insignificant. 

We may now proceed as above and eliminate all the 
irrelevant modes from the analysis by writing down the 
Gaussian Wigner function of the five interesting modes. 



The only difference is that this time (clcj) and {ciCj) are 
expressed in terms of the two time correlation functions 
of the OPO output. Also, the operators applied to the 
Wigner function to take conditioning into account are 
different because the annihilation detector model is used. 
The reader is referred to Refs. [l^ IISI for details. 

The resulting fidelity is shown as a function of e/7 in 
Fig.m where e is the nonlinear gain in the OPO, and as a 
function of the temporal distance between the condition- 
ing detection events in Fig. \S[ As in the pulsed case the 
fidelity decreases when the degree of squeezing increases. 
The fidelity also decreases when the temporal distance 
between the conditioning detection events increases from 
zero, but it is permissible to have a small time interval 
between the trigger detection events. We note that the 
curves represent a lower limit to the theoretically achiev- 
able fidelity since a better fidelity may be obtained for 
another choice of output state mode functions. 



V. CONCLUSION 

In conclusion we have analyzed a method to gener- 
ate path entangled NOON states from the output from 
two optical parametric oscillators. The method relies 
on the joint detection of photons in a number of trig- 
ger beams, and we presented a theoretical analysis of the 
role of detector efficiency and dead time for the fidelity 
of the states obtained and the success probability of the 
protocol. Our specific NOON state protocol applies for 
general photon numbers of the states, but in practice it 
is not realistic to go beyond the case of iV = 3, studied 
here. This is due to the reduction of the fidelity due to 
unwanted contributions from higher number states, when 
the OPO output power gets too high, combined with the 
severe reduction of the probability to obtain the number 
of conditioning detection events needed when the OPO 
output power is too low. The N — 3 NOON states, which 
can be produced at 90% fidelity at the rate of one state 
produced every 10 — 100 seconds, seem to be at the limit 
of realistic experiments of the proposed kind. Finally, we 
also determined the NOON state fidelity for continuous- 
wave fields, where the best NOON state occupies a pair 
of suitably selected temporal mode functions, and where 
we find high fidelities as long as the trigger events occur 
within a short time window compared to the lifetime of 
the OPO cavity field. 

We presented this analysis for the production of op- 
tical NOON states, but we note that recent theoreti- 
cal proposals and experiments with four wave mixing 
of matter waves |26|. e ngineered quadratic interactions 
among trapped ions [271 , and entanglement between field 
and atomic degrees of freedom ^28, 29] bring promise for 
similar conditional generation of atomic and interspecies 
atom-field NOON states. 

This work was supported by the European Integrated 
project SCALA. 
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